Baryon Stopping in Heavy-Ion Collisions at _E lab = 2—160 GeV/nucleon 
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It is argued that the experimentally observed baryon stopping may indicate (within the present 
experimental uncertainties) a non-monotonous behaviour as a function of the incident energy of 
colliding nuclei. This can be quantified by a midrapidity reduced curvature of the net-proton 
rapidity spectrum. The above non-monotonous behaviour reveals itself as a "zig-zag" irregularity 
in the excitation function of this curvature. The three-fluid dynamic calculations with a hadronic 
equation of state (EoS) fail to reproduce this irregularity. At the same time, the same calculations 
with an EoS involving a first-order phase transition into the quark-gluon phase do reproduce this 
"zig-zag" behaviour, however only qualitatively. 
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I. INTRODUCTION 

A degree of stopping of colliding nuclei is one of the 
basic characteristics of the collision dynamics, which de- 
termines a part of the incident energy of colliding nuclei 
deposited into produced fireball and hence into produc- 
tion of secondary particles. The deposited energy in its 
turn determines the nature (hadronic or quark-gluonic) 
of the produced fireball and thereby its subsequent evo- 
lution. Therefore, a proper reproduction of the baryon 
stopping is of prime importance for theoretical under- 
standing of the dynamics of the nuclear collisions. 

A direct measure of the baryon stopping is the net- 
baryon rapidity distribution. However, since experimen- 
tal information on neutrons is unavailable, we have to rely 
on proton data. Presently there exist extensive experi- 
mental data on proton (or net-proton) rapidity spectra 
at AGS [THI and SPS energies. These data were 
analyzed within various models II Oh 1811 The most exten- 
sive analysis has been done in [14l . |l7J • Since that time 
new data at SPS energies have appeared @-[|[. There- 
fore, it is appropriate to repeat this analysis of already 
extended data set. In the present Letter it is done within 
the framework of the model of the three-fluid dynamics 
(3FD) [13. 



II. ANALYSIS OF EXPERIMENTAL DATA 

Available data on the proton (at AGS energies) and 
net-proton (at SPS energies) rapidity distributions from 
central heavy-ion collisions are presented in Fig. [TJ Only 
the midrapidity region is displayed in Fig. [TJ since it is 
of prime interest in the present consideration. The data 
at 10 A GeV are repeated in the right panel of Fig. [JJ 
in order to keep the reference spectrum shape for the 



comparison. The data are plotted as functions of a "di- 
mensionless" rapidity (y - y cm )/y cm , where y cm is the 
ccnter-of-mass rapidity of colliding nuclei. In particular, 
this is the reason why the experimental distributions are 
multiplied y cm - This representation is chosen in order to 
make different distributions of approximately the same 
width and the same height. This is convenient for com- 
parison of shapes of these distributions. To make this 
comparison more quantitative, the data are fitted by a 
simple formula 
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where a, y s and w s are parameters of the fit. The form 
(UJ is a sum of two thermal sources shifted by ±y s from 
the midrapidity. The width w s of the sources can be 
interpreted as w s = (tcmperature)/(transverse mass), if 
we assume that collective velocities in the sources have 
no spread with respect to the source rapidities ±y s . The 
parameters of the two sources are identical (up to the 
sign of y s ) because we consider only collisions of identical 
nuclei. Results of these fits are demonstrated in Fig. 
[TJ Energy dependence of parameters y s and w s deduced 
from these fits revels no significant irregularities: they 
monotonously rise with the energy. 

The above fit has been done by the least-squares 
method. Data were fitted in the rapidity range \y — 
ycm\/ycm < 0.7. The choice of this range is dictated 
by the data. As a rule, the data are available in this 
rapidity range, sometimes the data range is even more 
narrow (40 A, 80 A GeV and new data at 158 A GeV 
We put the above restriction in order to treat different 
data in approximately the same rapidity range. Notice 
that the rapidity range should not be too wide in order 
to exclude contribution of cold spectators. 

We met problems with fitting the data at 80A GeV Q 
and the new data at 158 A GeV 0]. These data do not 
go beyond the side maxima in the rapidity distributions. 
The fit within such a narrow region results in the source 
rapidities y s very close (at 80A GeV) or even exceeding 
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FIG. 1: Rapidity spectra of protons (for AGS energies) and 
net-protons (p — p) (for SPS energies) from central collisions 
of Au+Au (AGS) and Pb+Pb (SPS). Experimental data are 
from collaborations E802 Q, E877 @|, E917 @], E866 [J, and 
NA49 .5 -9]. The percentage shows the fraction of the total 
reaction cross section, corresponding to experimental selec- 
tion of central events. Solid lines connecting points represent 
the two-source fits by Eq. ([1]). The dashed line is the fit to 
old data on Pb(158A GeV)+Pb 0, these data themselves are 
not displayed. 



(at 158 A GeV) y cm and a huge width w s . As a result, the 
normalization of the net-proton rapidity distributions, as 
calculated with fit (Ql, turns out to be 330 (at 80 A GeV) 
and 400 (at 158 A GeV), which are considerably larger 
than the total proton number in colliding nuclei (=164). 
To avoid this problem, we performed a biased fit of these 
data. An additional condition restricted the total nor- 
malization of distribution ([T]) to be less than the total 
proton number in colliding nuclei (=164). This biased fit 
is the reason why the curve fitted to the new data at 158 A 
GeV does not perfectly hit the experimental points. In 
particular, because of this problem we keep the old data 
at 158A GeV [f| in the analysis. We also use old data 
at 40^4 GeV, corresponding to centrality 7% Jj|, instead 
of recently published new data at higher (5%) centrality 
@, since the data at the neighboring energies of 20 A, 
30 A and 80 A GeV are known only at centrality 7% [|| . 
Similarity of conditions, at which the data were taken, 
prevents excitation functions, which are of prime interest 
here, from revealing artificial irregularities. 

Inspecting evolution of the spectrum shape with the in- 
cident energy rise, we observe an irregularity. Beginning 
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FIG. 2: Midrapidity reduced curvature of the (net)proton 
rapidity spectrum as a function of the center-of-mass energy 
of colliding nuclei as deduced from experimental data and 
predicted by 3FD calculations with hadronic EoS (hadr. EoS) 
[19l ] and a EoS involving a first-order phase transition into 
the quark-gluon phase (2-ph. EoS) [20t | ■ The thin dashed- 
dotted line demonstrates the effect of the 2-ph. EoS without 
changing the friction in the quark-gluon phase. 



from the lowest AGS energy to the top one the shape of 
the spectrum evolves from convex to slightly concave at 
10A GeV. However, at 20A GeV the shape again becomes 
distinctly convex. With the further energy rise the shape 
again transforms from the convex form to a highly con- 
cave one. In order to quantify this trend, we introduce 
a reduced curvature of the spectrum in the midrapidity 
defined as follows 
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This curvature is defined with respect to the "di- 
mensionless" rapidity (y — y C m)/ycm- The factor 
1/ (ycmdN I 'dy) y= y cm is introduced in order to get rid of 
overall normalization of the spectrum, i.e. of the a pa- 
rameter in terms of fit (JTJ). The second part of Eq. @ 
presents this curvature in terms of parameters of fit JT]). 

Values of the curvature C y deduced from fit (Q]) to ex- 
perimental data are displayed in Fig. [5J To evaluate er- 
rors of these deduced values, we estimated the errors pro- 
duced by the least-squares method, as well as performed 
fits in different the rapidity ranges: \y — y C m\/ycm < 0.5 
and \y — y cm \/y C m < 0.9, where it is appropriate, and also 
fits of the data at 80A GeV d| and the new data at 158 A 
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GeV 0] with different bias on the ovcrail normalization 
of the distributions: -/V prot < 208 (i.e., half of the net- 
nucleons can be participant protons) and N plot . < 128 
(which is the hydrodynamic normalization of the dis- 
tribution). The error bars present largest uncertainties 
among mentioned above. The lower point at s 1 / 2 = 17.3 
GeV corresponds to the new data at 158A GeV. Its up- 
per error, as well as that of 80^4-GcV point, results from 
the uncertainty of the normalization. The irregularity 
observed in Fig. Q]is distinctly seen here as a "zig-zag" 
irregularity in the energy dependence of C y . 

It is somewhat suspicious that the "zig-zag" irregu- 
larity happens at the border between the AGS and SPS 
energies. It could imply that this irregularity results from 
different ways of selecting central events in AGS and SPS 
experiments. However, there are indirect evidences of a 
physical (rather than methodical) nature of this irregu- 
larity. The difference between C y values in two different 
experiments at 10 A GeV can be taken as an estimate of 
the methodical uncertainty. The difference between C y 
values at 10A GeV and 20A GeV is two to three times 
larger than this methodical uncertainty. Moreover, we 
could expect that C y at 2QA GeV would be larger than 
that at 1QA GeV because the incident energy is higher 
and ccntrality selection at 2QA GeV is less restrictive 
(7%) than at 1QA GeV (5%). Contrary to these expec- 
tations the C y at 2QA GeV is smaller than that at 10A 
GeV. There should be a physical reason for that. Ex- 
citation functions of other quantities plj deduced from 
the same AGS and SPS data do not reveal any misht 
at the border between the AGS and SPS domains. The 
latter suggests that the AGS and SPS data were taken 
at similar physical conditions. However, new data taken 
at the same acceptance and the same centrality selection 
in this energy range are highly desirable to clarify this 
problem. Hopefully such data will come from new ac- 
celerators FAIR at GSI and NICA at Dubna, as well as 
from the low-energy-scan program at RHIC. 



III. THREE-FLUID MODEL SIMULATIONS 

Figure [5] also contains C y deduced from results of 3FD 
simulations with a hadronic equation of state (hadr. EoS) 
[l9| and a EoS involving a first-order phase transition 
into the quark-gluon phase (2-ph. EoS) [20| . To obtain 
y s and w s , the 3FD spectra were also fitted by the form 
([T]). For central (5%) Au+Au collisions at AGS energies 
we performed our calculations taking a fixed impact pa- 
rameter b = 2 fm; for the central (5%) Pb+Pb reaction 
at £i a b = 158A GeV, b = 2.4 fm which is the experimen- 
tal estimate for this centrality [22[ ; for other central (7%) 
Pb+Pb collisions at 20A-80A GeV, 6=3 fm. 

The 3FD model with the hadronic EoS reasonably re- 
produces a great body of experimental data in a wide 
energy range from AGS to SPS, see Ref. ^UMl. De- 
scription of the rapidity distributions with the hadronic 
EoS is reported in Refs. [UIHI- The reproduction of the 
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FIG. 3: Midrapidity value of the net-proton rapidity spec- 
trum as a function of the center-of-mass energy of colliding 
nuclei. Experimental data are confronted to predictions of 
the 3FD model with the hadronic EoS (hadr. EoS) pi] and 
the EoS with a first-order phase transition (2-ph. EoS) [20( |. 
The thin long-dashed line corresponds to the hadr. -EoS cal- 
culation without fragment production. 



distributions is quite good at the AGS energies and at the 
top SPS energies. At 40^4 GeV the description is still sat- 
isfactory. However, at 20A and 30A GeV the hadr. -EoS 
predictions completely disagree with the data, cf. fl8j . 
At 20A GeV instead of a bump at the midrapidity the 
hadronic scenario predicts a quite pronounced dip. The 
problems with the description of the low-energy SPS data 
are clearly seen from Fig. [2] and also from Fig. [3l where 
midrapidity values of the rapidity distributions (multi- 
plied by the center-of-mass rapidity) are presented. In 



Ref. [18j it was demonstrated that the problem with the 



low-energy SPS data can be solved by considerable soft- 
ening the hadronic EoS. This softening may indicate an 
onset of the phase transition into the quark-gluon phase. 
Notice that a maximum in y C m(dN/dy) cm at s 1 / 2 = 4.7 
GeV happens only because the light fragment production 
becomes negligible above this energy. The 3FD calcula- 
tion without coalescence (i.e. without the fragment pro- 
duction) reveals a monotonous decrease of y cm (dN/dy) cm 
beginning from s 1 / 2 = 2.7 GeV, i.e. from the lowest en- 
ergy considered here. 

The 3FD simulations have been also done with a EoS 
involving a first-order phase transition into the quark- 
gluon phase (2-ph. EoS) [20[. In 2-ph. EoS the Gibbs 
construction was used for the mixed phase. These calcu- 
lations well reproduce the AGS data up to the energy of 
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FIG. 4: The same as in Fig. [3] but in conventional repre- 
sentation (without multiplying by y C m) and in a wider energy 
range including RHIC data on Au+Au collisions at 5% cen- 
trality [2(|. Proton data 0, OSI are also displayed. 



6A GcV, where the purely hadronic scenario is realized. 
The data at the top SPS energy are also reproduced, 
which is achieved by a proper tune of the inter-fluid fric- 
tion in the quark-gluon phase. Quality of the reproduc- 
tion of above data is approximately the same as that 
with the hadronic EoS, as it is, e.g., seen from Figs. [5] 
and from Fig. |3] However, at top AGS and lower SPS 
energies (8A-80A GeV), where the mixed phase turns out 
to be really important, the 2-ph. EoS completely fails. 
The fact that the 2-ph.-EoS line perfectly hits 20A-4QA- 
GeV experimental points in Fig. [3] is just a coincidence, 
shapes of the distributions are completely wrong, as seen 
from Figs. [2] This failure cannot be cured by variations 
of neither the friction nor the freeze-out criterion. 

However, the C y curvature energy dependence in the 
first-order-transition scenario manifests qualitatively the 
same "zig-zag" irregularity (Fig. [5]), as that in the 
data fit, while the hadronic scenario produces purely 
monotonous behaviour. This "zig-zag" irregularity of 
the first-order-transition scenario is also reflected in the 
midrapidity values of the (net)proton rapidity spectrum 
(Fig. [3]). As for the experimental data, it is still diffi- 
cult to judge if the "zig-zag" anomaly in the midrapid- 
ity values is statistically significant. In the conventional 
representation of the data (Fig. [2J without multiplying 
by Verm the irregularity of the (dN/dy) cm data is hardly 
visible. However, the conventional representation clearly 
demonstrates the overall trend of the data: the midra- 
pidity net-proton yield gradually decreases with the in- 



cident energy, while the proton one stays approximately 
constant above the top SPS energy. Below the top SPS 
energy the proton and net-proton yields practically coin- 
cide. Model computations above the top SPS energy are 
at present not feasible because of high memory consump- 
tion required by the code (see discussion in Ref. [l?} )• 

All above discussion concerns only central nuclear colli- 
sions. Experimental data on midccntral collisions is much 
less complete. The model calculations for midcentral col- 
lisions (b w 6 fm) reveal the same quantitative behaviour 
of the excitation functions of C y and (dN/dy) cm both for 
hadr. EoS and 2-ph. EoS. 

The baryon stopping depends on a character of inter- 
actions (e.g., cross sections) of the matter constituents. 
If during the interpenetration stage of colliding nuclei a 
phase transformation 1 of the hadronic matter into quark- 
gluonic one happens, one can expect a change of the 
stopping power of the matter at this time span. This 
is a natural consequence of a change of the constituent 
content of the matter because hadron-hadron cross sec- 
tions differ from quark-quark, quark-gluon, etc. ones. 
This can naturally result in a non-monotonous behaviour 
of the shape of the (net)proton rapidity-spectrum at an 
incident energy, where onset of the phase transition oc- 
curs. Of course, the first-order transition does not hap- 
pen abruptly. Within the Gibbs construction the fraction 
of the quark-gluon phase is gradually increasing, as well 
as weights of the corresponding cross sections. Therefore, 
a non-monotonous behaviour will show up only if the dif- 
ference in cross sections in the hadronic and quark-gluon 
phases is large enough to override the above gradual in- 
crease of the fraction of the new phase. In fact, this is 
the case in the 3FD calculation with the phase transi- 
tion (2-ph. EoS). The friction in the quark-gluon phase 
was tuned to reproduce the data at the top SPS energy. 
Naturally, it does not continuously match the friction in 
the hadronic phase. In terms of parton-parton cross sec- 
tions, these cross sections in the quark-gluon phase turn 
out to be approximately twice as large as those in the 
hadronic phase 2 . In the quark-gluon phase these cross 
sections are compatible with those used, e.g., in a multi- 
phase transport model [27j and a parton cascade model 

Notice that the proton rapidity distribution at 158^4 
GeV is well described within the color-glass-condensate 
framework based on small-coupling QCD (29[. This 
mechanism drastically differs from that of hadronic stop- 
ping. Therefore, it is not surprising that the 3FD model 
requires very different (from the hadronic one) phe- 
nomenological friction at the 158A-GeV energy to repro- 
duce the data. 



1 The term "phase transition" is deliberately avoided, since it usu- 
ally implies thermal equilibrium. 

2 In the hadronic phase this parton cross section corresponds to 
the proton-proton one on the assumption of naive valence quark 
counting. 
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FIG. 5: Dynamical trajectories of the matter in the central 
box of the colliding nuclei (4fmx4fmX7 cm 4fm), where 7 cm is 
the Lorentz factor associated with the initial nuclear motion 
in the cm. frame, for central (b = 0) collisions of Au+Au 
at 4A and 10 A GeV energies and Pb+Pb at 20,4 GeV. The 
trajectories are plotted in terms of baryon density (ub) and 
the energy density minus ns multiplied by the nucleon mass 
(e — niNns)- Only expansion stages of the evolution are dis- 
played for two EoS's. Symbols on the trajectories indicate 
the time rate of the evolution: time span between marks is 1 
fm/c. 



However, if even the same friction is used in both 
phases, the calculated (with 2-ph. EoS) reduced cur- 
vature still reveals a "zig-zag" behaviour but with con- 
siderably smaller amplitude (see the thin dashed-dottcd 
line in Fig. [2]). This happens because the EoS in a gener- 
alized sense of this term, i.e. viewed as a partition of the 
total energy between kinetic and potential parts, also af- 
fects the stopping power. The friction is proportional to 
the relative velocity of the counter-streaming nuclei fl7j . 
Therefore, it is more efficient when the kinetic-energy 
part of the total energy is higher, i.e. when the EoS is 
softer. This effect of the softening was demonstrated in 
Ref. [lH . It was shown that application of a soft, but still 
hadronic EoS changes the rapidity distributions, making 
them closer to the data at low SPS energies. This is pre- 
cisely what the phase transition docs: it makes the EoS 
essentially softer in the mixed-phase region. The latter 
naturally results in a non-monotonous evolution of the 
proton rapidity spectra with the energy rise. 

Figure[5]demonstrates that the onset of the phase tran- 
sition in the calculations indeed happens at top-AGS- 



low-SPS energies, where the "zig-zag" irregularity takes 
place. Similarly to that it has been done in Ref. [30| . 
the figure displays dynamical trajectories of the matter 
in the central box placed around the origin r = (0,0,0) 
in the frame of equal velocities of colliding nuclei: \x\ < 2 
f m i I V | < 2 fm and \z\ < 2 fm, where 7 cm is Lorentz 
factor associated with the initial nuclear motion in the 
cm. frame. Initially, the colliding nuclei are placed sym- 
metrically with respect to the origin r = (0,0,0), z is 
the direction of the beam. The e-ns representation is 
chosen because these densities are dynamical quantities 
and, therefore, are suitable to compare calculations with 
different EoS's. Subtraction of the ra^qUB term is taken 
for the sake of suitable representation of the plot. Only 
expansion stages of the evolution are displayed, where 
the matter in the box is already thermally equilibrated. 
The size of the box was chosen to be large enough that 
the amount of matter in it can be representative to con- 
clude on the onset of the phase transition and to be 
small enough to consider the matter in it as a homo- 
geneous medium. Nevertheless, the matter in the box 
still amounts to a minor part of the total matter of col- 
liding nuclei. Therefore, only the minor part of the total 
matter undergoes the phase transition at 10A GeV en- 
ergy. As seen, the trajectories for two different EoS's are 
very similar at AGS energies and start to differ at SPS 
energies because of the effect of the phase transition. 



IV. CONCLUSIONS 

In conclusion, it is argued that the experimentally ob- 
served baryon stopping may indicate (within the present 
experimental uncertainties) a non-monotonous behaviour 
as a function of the incident energy of colliding nuclei. 
This reveals itself in a "zig-zag" irregularity in the exci- 
tation function of a midrapidity reduced curvature of the 
(net)proton rapidity spectrum. Notice that the energy 
location of this anomaly coincides with the previously 
observed anomalies for other hadron-production proper- 
ties at the low SPS energies [U [3l[ . The 3FD calculation 
with the hadronic EoS fails to reproduce this irregularity. 
At the same time, the same calculation with the EoS in- 
volving a first-order phase transition into the quark-gluon 
phase (within the Gibbs construction) [2(| reproduces 
this "zig-zag" behaviour, however only qualitatively. Pre- 
liminary simulations with the EoS of Ref. [32[ , also based 
on the first-order phase transition but within the Maxwell 
construction, show the same qualitative trend. It is ar- 
gued that the non-monotonous behaviour of the baryon 
stopping is a natural consequence of a phase transition. 
The question why these calculations do not qualitatively 
reproduce the "zig-zag" irregularity deserves special dis- 
cussion elsewhere. It is very probable that either the 
Gibbs and Maxwell constructions are inappropriate for 
the fast dynamics of the heavy-ion collisions [3^, [H| or 
the phase transition is not of the first order. 
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